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Abstract 

In recent years there has been growing interest in study of multi-antenna transmit designs for 
providing secure communication over the physical layer This paper considers the scenario of an intended 
multi-input single-output channel overheard by multiple multi-antenna eavesdroppers. Specifically, we 
address the transmit covariance optimization for secrecy-rate maximization (SRM) of that scenario. The 
challenge of this problem is that it is a nonconvex optimization problem. This paper shows that the SRM 
problem can actually be solved in a convex and tractable fashion, by recasting the SRM problem as a 
semidefinite program (SDP). The SRM problem we solve is under the premise of perfect channel state 
information (CSI). This paper also deals with the imperfect CSI case. We consider a worst-case robust 
SRM formulation under spherical CSI uncertainties, and we develop an optimal solution to it, again via 
SDR Moreover, our analysis reveals that transmit beamforming is generally the optimal transmit strategy 
for SRM of the considered scenario, for both the perfect and imperfect CSI cases. Simulation results 
are provided to illustrate the secrecy-rate performance gains of the proposed SDP solutions compared to 
some suboptimal transmit designs. 

Index terms— Physical-layer secrecy, secrecy capacity, transmit beamforming, semidefinite program. 

EDICS: MSP-CODR (MIMO precoder/decoder design), MSP-APPL (Applications of MIMO communi- 
cations and signal processing), SAM-BEAM (Applications of sensor and array multichannel processing) 



^ Copyright (c) 2011 IEEE. Personal use of this material is permitted. However, permission to use this material for any other 
purposes must be obtained from the IEEE by sending a request to pubs-permissions@ieee.org. This work is supported by a 
General Research Fund awarded by Research Grant Council, Hong Kong (Project No. CUHK415908). Part of this work appears 
in ICASSP 2010. 

^ Qiang Li is with Department of Electronic Engineering, The Chinese University of Hong Kong, Shatin, Hong Kong S.A.R., 
China. E-mail: qli@ee.cuhk.edu.hk. 

*Wing-Kin Ma is the corresponding author. Address: Department of Electronic Engineering, The Chinese University of Hong 
Kong, Shatin, Hong Kong S.A.R., China. E-mail: wkma@ieee.org. 



April 19, 2011 



DRAFT 



2 



I. Introduction 

Physical-layer secrecy is an information theoretic approach where we intend to provide a legitimate 
receiver with a reliable communication, and, at the same time, make sure that illegitimate receivers can 
retrieve almost nothing about the transmitted information from the signals they have intercepted. The 
study of this topic is meaningful and important, enabling us to understand the information rate limits 
when perfect secrecy is desired; i.e., the secrecy capacity or the maximum secrecy rate. Moreover, the 
physical-layer secrecy study provides us with vital implications on how physical-layer secret transmit 
schemes should be designed in practice. While the concepts of physical-layer secrecy can be found back 
in the 70's; e.g., the seminal works by Wyner [H, Leung- Yan-Cheong and Hellman El, and Csiszar and 
Korner ||3l, this topic has attracted much interest in recent years, in both information theory Bl- llTTI 
and signal processing |[T2l - |[23l . We can see at least two reasons for this. First, the rapid advances of 
wireless system architectures and applications, such as those for wireless networks, have given rise to 
new issues regarding information security. In particular, the open nature of the wireless medium means 
that signal interception may be easily conducted by eavesdroppers. Cryptographic encryption, the class of 
techniques commonly used to provide information security, is expected to be faced with more challenges; 
for instance, in key distribution and management |[24l . |[25l . Physical-layer secrecy suggests a physical- 
layer-based alternative to attacking the security problem, which is meaningful and may complement 
the network-layer-based cryptographic encryption techniques. Second, multi-input multi-output (MIMO) 
techniques provide physical-layer secrecy with new and exciting opportunities. Intuitively, if we know 
the channel state information (CSI) of the eavesdroppers to a certain extent (say, in scenarios where the 
eavesdroppers are users of the system, who attempt to access unauthorized services), then we may utilize 
the MIMO degree of freedom to weaken these eavesdroppers' receptions. Another idea that is possible 
only with MIMO is that of interfering eavesdroppers through artificially generated spatial noise; see ||T4l 
for the original work and |[T2l . |[T5l . ifTTll . ifTSl . |[23l for some recent developments. 

In fact, we have recently seen a growing body of physical-layer secrecy literatures that deal with various 
MIMO scenarios. For the scenario of an MIMO channel overheard by one multi-antenna eavesdropper, 
the secrecy capacity has been considered in Q-IU, |[26l . The MISO counterpart has also caught some 
attention; see HI, ifTTI . We should note that there is a difference between the above described MIMO and 
MISO scenarios, from a viewpoint of transmit optimization. Specifically, we are interested in transmit 
covariance designs for achieving the maximum secrecy rate. This secrecy-rate maximization (SRM) 
problem is nonconvex, and solving it under the general MIMO scenario is a challenging issue; see 



April 19, 2011 



DRAFT 



3 



mi, |[2ll for some recent endeavors. One exception where the SRM problem is found to be tractable is 
the MISO scenario — its optimal transmit design is shown to admit a closed-form solution, involving a 
generalized eigenvalue problem IH. 

This paper considers the scenario of an MISO channel overheard by multiple multi-antenna eavesdrop- 
pers. This scenario would happen, say, for example, in a downlink environment where the intended receiver 
uses one antenna for low operational costs and the eavesdroppers are users of the same network employing 
sophisticated multi-antenna hardware to improve their interceptions. The SRM problem for this particular 
scenario has been formulated in ifTOl . The problem is again nonconvex, and, to our knowledge, an efficient 
optimal transmit solution to it has not been available in general. We should mention that for the special 
case of one-antenna eavesdroppers, an optimal SRM design has been proposed in |[T6l : the idea there is to 
establish a relationship between the SRM problem and a cognitive radio (CR) design problem under the 
one-antenna eavesdroppers/secondary-users context. That secrecy-CR relationship result is applicable also 
to MIMO intended channel. For the multiple-antenna eavesdroppers case, the secrecy-CR relationship 
result does not apply and it is used as an approximation to the SRM problem |fT6l . This paper shows that 
the considered SRM problem can actually be solved in a convex and tractable fashion; specifically, we turn 
the SRM problem to a semidefinite program (SDP), a representative convex optimization problem whose 
globally optimal solution can be obtained efficiently by available algorithms |[27l . ||28]l . The proposed 
approach is indirect, and it may remind us of study of some kind of dual problem^!! in the context 
of multiuser downlink transmit beamforming |[29l . |[30l : also |[3ll . ll32l for multiuser downlink OFDM. 
In particular, we need to consider another secrecy-rate formulation, namely, a secrecy-rate constrained 
problem, and use the results established there to prove that the SRM problem has an SDP equivalent. 

Another key contribution of this work is that we extend our SRM results to the imperfect CSI 
case. The eavesdroppers' CSIs may be estimated depending on the application environment; e.g., when 
the eavesdroppers are also system users, however we may only have some inaccurate, possibly rough, 
knowledge about those CSIs. The CSI of the intended receiver may also be subject to some uncertainties; 
e.g., those caused by channel estimation errors and/or quantization errors, though such an issue should 
be less severe than that for the eavesdroppers. We should note that physical-layer secrecy with imperfect 
CSI or no CSI has started to catch attention very recently; see, e.g., H, IBl, HH, lH, ||33l. Here, 
we consider a robust SRM formulation of our considered scenario, where we employ a spherical CSI 

'The dual problem relationship we refer to is that of the quality-of-service (QoS) constrained power minimization problem 
and the power constrained QoS maximization problem 1291 . 
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uncertainty model and the uncertainties are handled in a worst-case sense. The proposed robust SRM 
formulation is a conservative design — it guarantees perfect secrecy for any admissible CSI uncertainties, 
including the worst. The robust SRM formulation exhibits a more complex structure than its non-robust 
counterpart. Fortunately, it is shown that the robust SRM problem can be equivalently represented, and 
solved, by an SDP; in addition to the approach used in non-robust SRM, the idea is to employ linear 
matrix inequality characterization to handle the imperfect-CSI induced constraints in a tractable way. 

In establishing our optimal SDP solutions to the SRM problems, we obtain a physical result that agrees 
well with intuitive expectation — transmit beamforming is generally the SRM optimal transmit strategy of 
the considered scenario, for both the perfect and imperfect CSI cases. While this result is already known 
for the one multi-antenna eavesdropper case H and the multiple one-antenna eavesdroppers case |[T6l 
(both with perfect CSI), the result here is more general. 

This paper is organized as follows. A background review and problem statement is given in Section lU 
Section |lll] considers the SRM problem for the MISO, multi-eavesdropper scenario with perfect CSI, 
and establishes an SDP solution to it. Section |IV] extends the SRM results to the imperfect CSI case. 
Simulation results comparing the proposed SRM solutions and some other suboptimal secrecy transmit 
designs are illustrated in SectiorfVl Section |Vl] gives the conclusion and discussion. 

Our notations are as follows. We will use boldface capital letters to denote matrices, boldface lower case 
letters to denote vectors; A'^, A^, A^, Tr(A) and det(A) represent transpose, Hermitian (conjugate) 
transpose, pseudo inverse, trace and determinant of a matrix A; I denotes an identity matrix; || • \\f 
represents the Frobenius norm of a matrix; rank(A) is the rank of matrix A; A ^ (A ^ 0) means A 
is a Hermitian positive semidefinite (definite) matrix; denotes the set of all nonnegative real numbers; 

denotes the set of all A^-by-A^ Hermitian matrices; denotes the set of all A^-by-A^ Hermitian 
positive semidefinite matrices; C^^^^ represents an M-by-A^ dimensional complex matrix set; A ig) B 
denotes the Kronecker product of A and B; vec(A) denotes the vectorization of matrix A by stacking 
its columns; E{-} is the expectation operator 

II. Background and Problem Statement 

In this section we review some basic concepts of physical-layer secrecy with an emphasis on the MISO 
scenarios, and provide the problem statement of our interested scenario. 
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A. One Eavesdropper Case: A Review 

The endeavor of this subsection is to give some intuitive insights into physical-layer secrecy, by 
reviewing the scenario of a MISO channel overheard by one eavesdropper (also known as the MISOME 
wiretap channel in the literature pP]). Fig. [Ha) shows the scenario. In this configuration, the transmitter 
is equipped with multiple antennas, intending to deliver information to the legitimate receiver which uses 
one antenna. One eavesdropper is present, and it uses multiple antennas to intercept the transmission. 
Following the convention in the secrecy capacity literature, we call the transmitter, the legitimate receiver, 
and the eavesdropper as Alice, Bob, and Eve, respectively. The signal models for the Alice-to-Bob and 
Ahce-to-Eve links are, respectively, 

= h^x(t)+n(t), (la) 
ye(t) = G^x(t) + v(t). (lb) 

Here, x(t) G C^* is the transmit vector by Alice, with Nt being the number of transmit antennas; h G C^' 
represents the MISO Alice-to-Bob channel; G G ([^Nty.N^ represents the MIMO Alice-to-Eve channel, 
with A'^e being the number of receive antennas at Eve; and n{t) G C, v(t) G C^" are zero-mean additive 
white Gaussian noises. Without loss of generality, we assume unit variance of all the noise terms; i.e., 
E{Kt)|2} = l, E{v(t)v^(t)}=I. 




(a) (b) 
Fig. 1. System model, (a) One multi-antenna eavesdropper; (b) Multiple multi-antenna eavesdroppers. 

The problem of interest here is to prevent Eve from retrieving useful information via appropriate 
transmit designs, and this can be addressed by using the notion of physical-layer secrecy |[34l . To describe 
this, let us denote the transmit covariance by 

W = E{x(t)x^(t)}. 
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According to ||4l, the transmit covariance design that provides the maximum secrecy rate for the system 
model © is given by 

R^(P) = max log(l + h^Wh) - log detfl + G^WG) 

w (2) 

s.t. W ^ 0, Tr(W) < P, 
where P is a given average transmit power limit, and R*{P) is defined to be the optimal secrecy 
rate (in bps/Hz) for a given P. As seen in Q, the problem is to maximize the mutual information 
difference between the Alice-to-Bob and Alice-to-Eve channels. The subsequent optimized rate R*{P) 
is achievable — from an information theoretic perspective, there exist codes such that Bob can obtain a 
perfectly secure message from Alice at a rate of R*{P) bps/Hz, while Eve can retrieve almost nothing 
about the message |[3]. It has also been shown that Q is the secrecy capacity for the one-Eve model 
in ©. 

The secrecy-rate maximization problem in Q has a closed-form solution. Let q be the unit-norm 
principal generalized eigenvector of (I+Phh^, I+PGG^). The optimal transmit design of Q, denoted 
herein by W*, can be shown to be lU 

^ f Pqq^, /(Pqq^) > 
W* = <^ (3) 

0, /(Pqq^)<0 

V 

where /(W) = log(l + h^Wh) — logdet(I + G^WG) denotes the secrecy rate of a given W. The 
structure of W* shown above reveals two physical results: First, the optimal transmit strategy for the 
case of R*{P) > (positive secrecy rate) is to employ a rank-one transmit structure; i.e., transmit 
beamforming. Second, the optimal transmit design for the (trivial) case of R*{P) = is to shut down 
the transmitter. 

B. Multi-Eavesdropper Case 

This paper focuses on the scenario of MISO Alice-to-Bob link and multiple MIMO Alice-to-Eves 
links; such a configuration is depicted in Fig. [Hb). For this scenario the system model is modified to 

= h^x(t)+n(i), (4a) 
ye,fc(t) = Gfx(i)+vfc(t), k = l,...,K, (4b) 

where the Alice-to-Bob channel model in (|4al ) is identical to that in ([Tat ; ye,fe(i) is the received signal at 
the kth Eve; G^ E c^^txAf^.i (-j^e MIMO channel from Alice to the kth Eve; N^^k is the number of receive 
antennas of the kth Eve; Vk{t) is additive white Gaussian noise with zero mean and E{vfc(t)vfc(f)^} = I; 
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K is the number of Eves. An achievable secrecy rate for the multiple -Eve model (01) has been derived 
in ifTOl and is given by 

R^iP) = max mill /fc(W) 

W k=l,...,K (-5-) 

s.t. W ^ 0, Tr(W) < P 

where 

fk(W) = log(l + h^Wh) - log det(I + Gf WGfe) (6) 

is the secrecy rate function corresponding to the fcth Eve. The goal of ^ is to maximize the worst 
secrecy rate (or mutual information difference) among all the Eves. 

The secrecy-rate maximization (SRM) problem ^ presents a challenge from the standpoint of transmit 
covariance optimization. Problem (l5]l is a nonconvex problem, due to the nonconcave, Eve-induced terms 
- log det(I + G^WGk). While this nonconvex SRM problem admits a closed-form solution in the one- 
Eve case (as mentioned above), it is not known if the SRM problem can also be solved in closed form 
for the general multi-Eve case. 

We should mention a transmit design that has a closed-form solution but is generally suboptimal in 
the SRM context, namely, projected maximum-ratio transmission (projected-MRT). The idea is to apply 
nulling on all Eves; i.e., to enforce G|^W = for all k. Let G = [ Gi, . . . , Gk ] be the aggregate 
channel matrix of all Eves, and IIq = I — G{G^ G)^G^ be the orthogonal complement projector of 
G. The transmit beamformer weight in projected-MRT is 



y P 

l|n^h||2 

or, in other words, we choose W = ww^ . The advantage of projected-MRT lies in its simplicity, but 
its downside is that we may not have the degree of freedom to perform nulling when the total number 
of antennas of Eves, Xlfe^i ^e,k, reaches or exceeds the number of transmit antennas Nt. 
In what follows, we will describe our approach to solving the SRM problem Q. 

III. An SDP Approach to the Secrecy Rate Maximization Problem 

The focus of this section is on solving Problem ((Sjl, the secrecy-rate maximization of an MISO channel 
eavesdropped by multiple multi-antenna eavesdroppers. We will show how the SRM problem ([5]), which 
is nonconvex, can actually be solved by an equivalent SDP problem that is convex and tractable. 

The idea behind the proposed SDP solution is to consider some form of convex approximation to the 
SRM and then to prove that that approximate SRM problem is indeed tight. The latter part is not 
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Straightforward. We will need to study a variation of SRM, given as follows: 

P*{R) = min Tr(W) 
Who 

s.t. min /fc(W) > R. 

k=l,...,K 

The design in dV) seeks to satisfy a minimum secrecy rate specification R (which is given) and to 
minimize the average power. The reason for considering the secrecy-rate constrained (SRC) problem ^ 
is that it is relatively easier to analyze than the SRM (H). And, as a side advantage, the SRC design 
itself is interesting and practically meaningful. We will develop some key results for Problem (|7]), e.g., 
that (|7]) has an SDP equivalent. We will then establish a link between the SRC problem ^ and the SRM 
problem (jSj, thereby allowing us to use the proven results in SRC to show that SRM can be exactly 
solved by an SDP. 

The first and second subsections describe our developments for the SRC and SRM problems, respec- 
tively. 

A. The Secrecy-Rate Constrained Problem 

To study the SRC problem d?]), let us express ([7]) in a more explicit form: 

P*{R) = min Tr(W) 

R det(I + G^WGfc) 

s.t. 2-^> max \ , ^' . 

- k=i,...,K 1 + h^Wh 

Problem ([8]l is nonconvex, due to the determinant functions in the constraint of ([8]l. Consider the following 
lemma which we will use to provide a convex approximation to 

Lemma 1 Let A y 0. It holds true that 

det(I + A) > 1 + Tr(A), (9) 

and that the equality in ^ holds if and only if rank( A) < 1. 

Proof: Let r = rank(A). The case of r = is trivial. For r > 1, let Ai > A2 > • • • > > denote 
the non-zero eigenvalues of A. We have that 

det(I + A) = nLi(l + Ai) = 1 + EI=i A. + Ei^fc AiAfc + . . . 

> 1 + ELi A, = l + Tr(A) 
and it can be seen that the equality above holds if and only if r = 1. ■ 
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Applying Lemma 1 to Problem ([8]l, we obtain a relaxation of ([D as follows: 

1 + Tr(Gl^WGfc) (10) 



P'{R) > P:,iaAR) = min Tr(W) 
w >-o 



(11) 



s.t. 2 ^ > max ^ , 

~ k=i,...,K l + h-f^Wh 

Problem (fTOl ) is a convex problem. Specifically, (fTOb can be formulated as an SDP 

n*e^a.(^)=mm Tr(W) 

s.t. W >i 0, 

1 + Tr(hh^W) > 2^(1 + Tr(GfcGf W)), 
k = l,...,K, 

whose globally optimal solution can be efficiently found by available solvers ll27l . |[28l . While our 
original motivation is to approximate the SRC problem ([8]l by the convex relaxation ( fTOl) . Lemma 1 
provides an important hint regarding the tightness of the relaxation: The relaxation (ITOl ) is tight (which 
means P*{R) = P*eiax(^)^ when the optimal solution of (ITOl ) is of rank one. We prove the following 
key result: 

Proposition 1 Consider the relaxed SRC problem dlOl ) for the case where the secrecy-rate specification 
R is positive ^ or simply R > 0. Also, suppose that Problem (1101 ) is feasible. Then, the optimal solution 
of (1101 ) must be of rank one and unique. 

The proof, which can be found in Appendix A, is based on examination of the Karush-Kuhn-Tucker 
(KKT) conditions of the SDP (EJ. 

Using Proposition 1 and Lemma 1, we reach the conclusion that 

Corollary 1 Consider the SRC problem ^for the case R > 0, and suppose that Problem ([8]) is feasible. 
The relaxed SRC problem (1101 ) [ or the SDP (II II ) ] exactly solves the SRC problem, in the sense that the 
optimal solution of (1101 ) is also that of ([8]), and vice versa. Moreover, the optimal SRC solution is unique 
and of rank one. 

Proof: Let W be the optimal solution of (ITOl ). which is unique and of rank one by Proposition 1. By 
Lemma 1, W is feasible to ([8]), which implies that Preiaxi^) ~ Tr(W) > P*{R). Since we also have 
P*{R) > P*^iaAR) [cf., Eqn. we conclude that P*^iaAR) = P*{R); i.e., W is optimal to ^. 

^The case of _R = is trivial since W = can easily be verified to be the corresponding optimal solution of l llOt . 
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On the other hand, let W* be an optimal solution of ([D. Owing to Lemma 1, W* is feasible to (ITOl ). 
The result -P*e,„^(i?) = = Tr(W*) further implies that W* is optimal to (HOjl, too. By Proposition 

1, W* has to be unique and of rank one. ■ 

We have developed an SDP solution to the SRC problem ([8]l. Moreover, the unique rank-one nature 
of the SRC solution, as indicated in Corollary 1, suggests that transmit beamforming is generally the 
optimal transmit strategy for the SRC problem ([8]l. 

B. The Secrecy-Rate Maximization Problem 

We now turn our attention to the SRM problem ([5]). For convenience, we rewrite Problem ^ as 

7* P = min max \ , (12) 

' ^ ^ Who k=i,...,K l + h-f^Wh 

Tr(W)<P 

where < 7*(P) < 1 is related to the optimal secrecy rate through the relation R*{P) = log(l/7*(P)). 
Following the same spirit as in the preceding subsection, we apply Lemma 1 to (fT2l ) to obtain a relaxation: 

7* P > ireiaxiP) = ^ax ^ \ . (13) 

J \ J - lrelax\ I fc=l,...,K 1 + h^Wh 

Tr(W)<P 

The relaxation above is a quasi-convex problem, whose globally optimal solution can be searched by 
general techniques for quasi-convex optimization (e.g., bisection |[35|). We however will propose a more 
efficient method of solving ( fT3] ); namely, via SDP, after we study the tightness of the relaxation in ( fT3] ). 

Our key question here is whether Problem (fT3] ) yields a rank-one solution: If it does, then the 
relaxation ([T3] ) is tight (7*(P) = Ireiaxi-^)^ Lemma 1. Our insight to this is as follows: Suppose that 
^reiax(P) G^) has been computed. We consider the following relaxed SRC problem 

min Tr(W) 
w>rO 

. . s l + Tr(G^WGfe) ^^^> 

irelax\ k=l,...,K 1 + h-^Wh 

where we fix the secrecy rate specification at i? = log(l/7*g^^^(P)). By Proposition 1, we know that 
Problem ([T4b has a unique rank-one solution (in the nontrivial case). If the optimal solution of ([T4l ) can 
be proven to be an optimal solution of ([T3] ) as well, then we will be able to infer immediately that the 
relaxation ( [T3] | is tight. With this idea in mind, we give a formal proof in Appendix B and show that 
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Theorem 1 Consider the SRM problem for the case of < 7*(P) < 1. The relaxed SRM 
problem (1131 ) exactly solves the SRM problem, in the sense that the optimal solution of ([13b is also 
that of (fT2l) . and vice versa. Moreover, the optimal SRM solution is unique and of rank one^ 



The remaining issue is to solve the equivalent SRM problem (fT3l) . Problem (fT3l ) can be solved by 
using a bisection search methodology commonly used in quasi-convex optimization. However, such a 
method would require solving a sequence of (usually many) SDPs. Here we propose a more efficient 
alternative where we exploit the problem structures and recast (fT3] ) as an SDP. The idea is to use the 
Chames-Cooper transformation |[36l . |[37l . Let us consider the following transformation of the transmit 
covariance: 

W = Z/^ 



for some Z ^ 0, ^ > 0. The change of variables above enables us to re-express ( [13]) as 

^ + Tr(GfcGfZ) 
mm max — — ^ ,^ , (15a) 
z,5 k=i,...,K ^ + Tr(hh^Z) ^ ^ 

s.t. Tr(Z) < ^P, (15b) 

Z >i 0, e > 0, (15c) 
which may further be reformulated as an SDP 

min r (16a) 

S.t. C + T¥(GfcGf Z) < r, k = l,...,K, (16b) 

C + Tr(hh^Z) = 1, (16c) 

Tr(Z) < ^P, (16d) 

Z >r 0, ^ > 0. (16e) 

Note that (I16al )- (ll6bl ) is a consequence of the standard epigraph reformulation (see the literature, e.g., |[35l ). 
and the constraint (I16cl ) is additionally introduced to fix the denominator of the objective function in (I15al ) 
which is without loss of generality. The relatively more crucial part with the reformulation lies in replacing 
^ > in (I15cl ) by ^ > in (I16el) . We show that this will not cause a problem: 

^The case of ■J*{P) = 1 is trivial because that corresponds to zero secrecy rate and the respective SRM solution is simply 
W = 0. 
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Proposition 2 The SDP (1161 ) is equivalent to Problem (115b . The former is also equivalent to the SRM 
problem (|12|) through the relation W = Z/^. 

Proof: The remaining, nontrivial part is when ^ = in ( |16e| ). Suppose that this is true. Then, by ( |16d| ) 
and Z ^ 0, we must have Z = 0. The constraint (I16cl ) is then violated as a consequence. Thus, a feasible 
point of ([T6l ) must not have C = 0- This prove that Problem ( [T6l ) is equivalent to Problem ( fTSl ). The 
solution equivalence of Problems ( fTSl ) and ([T2l ) through W = Z/^ follows from the discussion above 
and Theorem 1. ■ 



Concluding, we have shown that the SRM problem ([T2l ) can essentially be solved by using SDP; see 
Theorem 1 and Proposition 2. Moreover, Theorem 1 reveals that the SRM solution must be of rank one, 
which implies that transmit beamforming is generally the optimal transmit strategy for the SRM problem. 

IV. Secrecy-Rate Optimization with Imperfect CSI 

The MISO secrecy problems solved in the previous sections have been based on a premise that the 
CSIs of Bob and Eves are perfectly known to Alice. In this section, we extend our results to the imperfect 
CSI case. We will consider robust MISO secrecy-rate formulations that cater for CSI uncertainties in the 
worst-case sense. The proposed robust secrecy-rate problems are more complex in structures and more 
challenging to solve than their perfect-CSI counterparts, but we will show that these problems can still be 
turned to SDPs. The robust formulations will be presented in the first subsection, while the SDP solutions 
to the robust formulations will be described in the second and third subsections. 

A. Robust Secrecy-Rate Problem Formulations 

Our model assumption for imperfect CSI is based on the deterministic model |[38l - H0l . We model the 
Alice-to-Bob channel and the Alice-to-Eve channels respectively by 

h = h + Ah, 

Gfc = Gfc + AGfc, k = l,...,K. 

Here, h and Gfc are the channel means and they are known to Alice; Ah and AG^ represent the 
channel uncertainties. These uncertainties are assumed to be deterministic unknowns with bounds on 
their magnitudes: 

||Ah||2 = ||h - h||2 < efe, 

||^Gfc||_F = ||Gfc - Gfclli? < Ee^fc, k = l,...,K, 
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for some Sb, ee,i, • • • , £e,K > 0. 

The proposed robust SRM formulation is as follows: 

R^iP) = max V'(W) (17) 

Tr(W)<P 

where i/'(W) is the worst-case secrecy rate function, defined by 

V'(W) = ^min^^fc(W), (18) 
V'fc ( W) = min log ( 1 + Wh) - log det (I + Gf WG^ ) , (19) 

Bb = {h£C^* \ ||h-h||2 <efe }, (20) 
B,,k = {Gk€ C^'^^-'= I \\Gk - Gk\\F < ee,fc }■ (21) 

The function ipk{-) represents the worst secrecy-rate function among all channel possibilities, for the 
kth Eve. The resulting design ([TT] ) is a conservative one; from its formulation we see that the secrecy 
rate will be guaranteed to be no less than the worst-case optimum R*{P) for any channel possibilities 
(described by Bb, Be^k)- 

Based on the same philosophy as in the last section, we can also formulate a worst-case robust SRC 
design: 

P^iR) = min Tr(W) 

w^o (^22) 

s.t. V(W) > R, 

i.e., minimizing the transmit power while satisfying a secrecy rate specification R in the worst-case sense. 
This robust SRC problem is interesting in its own rights, and, like the perfect-CSI solution established in 
the last section, solving the robust SRC problem will provide us with a crucial key to solving the robust 
SRM problem. 

B. The Robust Secrecy-Rate Constrained Problem 

We consider the robust SRC (R-SRC) problem (l22l ) in this subsection, where we aim to develop an 
SDP solution to R-SRC optimization. From (fT8])-(l22l). Problem (l22l ) can be expressed as 



P*iR) = min Tr(W) 



max det(I + GfWGfc) 



s.t. 2~''>^^ J, , (23) 

min 1 + h^Wh 

k = l,...,K. 
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Our approach to solving (l23l) is somehow reminiscent of its non-robust counterpart in Section IIII-AI We 
use Lemma 1 to relax the constraints in (l23l) . and obtain the following relaxed problem: 



P^R) > Pr\i,AR) = mm Tr(W) 

max l + Tr(GfWGfc) 

s.t. 2-''>^^ n , (24) 

mill 1 + h^Wh 

k = l,...,K. 

Our goals are then to turn (l24l ) to a tractable convex problem, and to show that the relaxation in (l24b is 
tight by proving the rank-one solution structure of (l24l ). 

In fact, the relaxed R-SRC problem (l24b can be recast as an SDR To do so, the first step is to 
reformulate (l24l) as 

min Tr(W) (25a) 

WbO,6» 

S.t. min 1 + h^Wh > 9, (25b) 
2-^e> max l + Tr(GfWGfc), (25c) 
k = l,...,K, 

where we have added a slack variable 9 to decouple the fractional functions in the original constraints 
in (l24l ). Problem ( |25] ) is already a convex problem in principle, but with semi-infinite constraints as seen 
in (|25b| ) and ( |25c| ). To make the problem more tractable to solve and analyze, the second step is to 
turn (I25bl ) and (I25cl ) to linear matrix inequalities (LMIs), using the 5-procedure: 

Lemma 2 ( 5-procedure Il35l ) Let 

/fe(x) = x^AfcX + 2Re{h^x} + 

for k = 1,2, where G H", G C", G M. T/ie implication /i(x) < =^ /2(x) < /loZt/s if and 
only if there exists /i > such that 

Ai bi A^ b9 

^0, 





Ai 


bi 




A2 b2 












bf 


Cl 




bf C2 



provided that there exists a point x smc/j that /i(x) < 0. 

To apply the 5-procedure to the constraint (I25bl ). we substitute the representation h = h+Ah into (I25bb 
and re-express (I25bl ) as: 

Ah^Ah < ^ Ah^WAh + 2Re{h^WAh} + h^Wh + 1 - 6* > 0. (26) 

April 19, 2011 DRAFT 



15 



Using the 5-procedure, we transform (1261 ) to an LMI 



Wh 

-Xbel -e + h^Wh + 1 



^0, 



(27) 



for some A5 > 0. Similarly, ( |25c| ) is equivalent to the following implication: 

Agj^Agk < elk =^ Agj^Wk^gk + 2Re{gj^WkAgk} + Wfc^fc + 1 - 2~^e < 0, (28) 

where Wfc = Ia^^ ^ (8) W and gk = vec(Gfc). By the 5-procedure, the above implication can be re- 
expressed as the following LMI: 

Ae,fcl7Ve,fcAft - Wfc -WkSlk 

for some Xe^k > 0, k = 1, . . . , K. Substituting (1271 ) and ( |29l ) back into (l25l ). we obtain the following 
SDP 



Te,fc(W,Ae,fc,^) 



h 0, 



(29) 



^'^eiaxl^) = min Tr(W) 



S.t. T;,(W,Afe,0) ^ 0, 

Te,fc(W,Ae,fc,0) ^0, k = l,...,K, 

W ^ 0, Aft > 0, Ae fe > 0, k = l,...,K. 



(30a) 

(30b) 
(30c) 
(30d) 



where Ag = [ Ae,i, . . . , Xe,K f- The SDP as an equivalent form of the relaxed R-SRC problem, 
can be solved conveniently by available SDP solvers. In addition to this, the problem structures of (l30l ) 
enable us to analyze the solution optimality of the relaxed R-SRC problem via the KKT conditions: 



Proposition 3 Consider the relaxed R-SRC problem (1301 ) for the case of R > 0, Also, suppose that 
Problem (1301 ) is feasible. Then, the optimal solution of (1301 ) must be of rank one and unique. 



It is interesting to note that the end results of Proposition 3 are exactly the same as those of its non- 
robust counterpart. Proposition 1. However, to prove Proposition 3 is considerably more complicated, 
owing to the more complex LMIs in (I30bl ) and (I30cb . The proof of Proposition 3 is given in Appendix 
C. 

We complete this subsection by using Proposition 3 to verify that the relaxed R-SRC problem (l30l ) is 
tight: 



Corollary 2 Consider the R-SRC problem (1231 ) for the case of R > 0, and suppose that problem (1231 ) 
is feasible. The relaxed R-SRC problem (1301 ) exactly solves the R-SRC problem, in the sense that the 
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optimal solutions of (1231 ) and (1301 ) are equivalent to one another. Moreover, the optimal R-SRC solution 
is unique and of rank one. 

Proof: The proof is essentially identical to that of Corollary 1, and thus is omitted for brevity. ■ 



C. The Robust Secrecy-Rate Maximization Problem 

With the R-SRC results established in the previous subsection, we are now ready to develop an SDP 
solution to the robust SRM (R-SRM) problem ( [TT] ). Problem ([TTl ) can be expressed as 

max det(I + Gl^WGfc) 

7 [F] = mm max (31) 

w>rO, k=i,...,K mill 1 + h^^ Wh 

Tr(W)<P heBb 

where < 7*(-P) < 1. Applying Lemma 1 to (|3TI ) yields the following relaxed R-SRM problem 

max l + Tr(GfWGfc) 

l\P)>lUAP)= min max ^ -\ , ^.H^^ • (^2) 

w^o, k=i,...,K mm 1 + h Wh 
Tr(W)<P heBi, 

We again investigate two issues: the tightness of the relaxed R-SRM problem, and the possibility of 
converting (|32l ) to a convex problem. 

The relaxed R-SRM problem is tight. We show the following: 

Theorem 2 Consider the R-SRM problem for the case of < j*{P) < 1. The relaxed R-SRM 
problem (1321 ) exactly solves the R-SRM problem, in the sense that the optimal solutions of (1311 ) and (1321 ) 
are equivalent to one another. Moreover, the optimal R-SRM solution is unique and of rank one. 

Proof: The proof of Theorem 2 is essentially the same as that of its non-robust counterpart. Theorem 1 , 
and here we provide only the outline. The idea is to consider the relaxed R-SRC problem 

min Tr(W) 

max l + Tr(GfWGfc) 

s.t. ll^^P) ^^^n^^ , (33) 

mm 1 + h W h 

k = l,...,K. 

It is shown that the optimal solutions of (l32l ) and (1331 ) are equivalent to one another, by following the 
same procedure as in Appendix B. The solution equivalence of (l32l ) and (1331 ). together with Proposition 
3, enable us to deduce that (l32l ) has a unique rank-one solution. In the same spirit as the proof of 
Corollary 1, the unique rank-one solution characteristic of ( [32] ) further imphes that 7*(-P) = l^eiax^P) 
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(recall Lemma 1), and that the solution of (l32l ) has to be the solution of (|3T]) . and vice versa. ■ 

The relaxed (actually, equivalent) R-SRM problem (l32l ) can be transformed to an SDR By employing 
the Charnes-Cooper transformation and 5-procedure, and by some careful manipulations, we show that 

Proposition 4 Problem (l32l ) is equivalent to the following SDP 

min r (34a) 

s.t. Mb(Z,Ab,e) >=0, (34b) 
Me,fc(Z, Ae,fc, r) ^ 0, k = l,...,K, (34c) 
TV(Z) < (34d) 
Z ^ 0, ^ > 0, Afo > 0, Ae,fe > 0, A; = 1, . . . , E:, (34e) 

where Xe = [ Ae,i, . . . , Xe,K 

Afoliv, + Z Zh 

h^Z h^Zh + e - Abe2 - 1 

-9k ^k -K,ke%k + T - Qk^kQk 
and Zk = lAfe.fc ® Z. Specifically, Problems (l32ll a«J (l34l) are equivalent through the relation W = Z/^. 

We delegate the proof of Proposition |4] to Appendix |Dl since the key ideas behind the proof, the 
Chames-Cooper transformation and 5-procedure, have been demonstrated in the preceding development. 

In summary, we have shown that the nonconvex R-SRM problem (|3T] ) can be equivalently solved by 
solving the convex SDP (l34l ): see Theorem 2 and Proposition 4. 

Recall that in the perfect CSI scenario studied in the previous section, transmit beamforming is shown 
to be the optimal transmit strategy for the SRC and SRM designs in general. This physical result remains 
valid for the worst-case robust SRC and SRM designs, as indicated in Corollary 2 and Theorem 2. 

V. Simulation Results 

We provide simulation results to illustrate the secrecy-rate performance gains of the proposed SDP 
solutions compared to some other existing methods. We will first consider the perfect CSI case in the 
first subsection, and then the imperfect CSI case in the second subsection. 



M6(Z,Afe,e) = 



Me,fc(Z,Ae,fc,e,r) = 
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A. The Perfect CSI Case 

The results to be presented in this subsection are based on the following simulation settings, unless 
specified: At Alice, the number of transmit antennas is Nt = 10, and the average transmit power limit 
is P = 3dB. Three Eves are present (K = 3), and their numbers of receive antennas are A'^e i = . . . = 
^6,^ = 3. Perfect CSIs are assumed. At each trial of the simulations, Bob's channel h is randomly 
generated following an independent and identically distributed (i.i.d.) complex Gaussian distribution with 
zero mean and unit variance. Similarly, each Eve's channel is randomly generated following an i.i.d. 
complex Gaussian distribution with zero mean and variance p^. We fix pg = 1, if not mentioned. The 
simulation results were obtained based on an average of 1000 independent trials. 

The following transmit designs are tested in our simulations: the proposed SDP solution ( fT6l ) to solving 
the SRM problem ([5]), the projected-MRT method described in Section ITl-B[ and a simple method called 
plain-MRT here. In plain-MRT we choose Bob's channel direction, h/||h||, as the transmit weight. 
Specifically, plain-MRT sets W = (P/||h|p)hh^ if the corresponding secrecy rate is positive; and 
W = otherwise. Plain-MRT is a suboptimal, arguably weak, method since it ignores the presence of 
Eves. It is nonetheless interesting to examine its secrecy-rate performance relative to SDP. 




Fig. 2. Secrecy rates of the various metiiods versus (a) the number of Eves, and (b) the average channel strength of Eves. 

1 ) Secrecy rates versus the number of Eves: Fig. |2fa) shows the secrecy rate behaviors of the various 
methods when we increase the number of Eves K. We can see that the proposed SDP method yields 
better performance than the two other methods over the whole range of K tested. The secrecy rate of 
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projected-MRT is able to approach that of SDP for K <2, but becomes zero for K > the latter case is 
when the degree of freedom of all Eves combined, XlfcLi ^e,k = 3i^, is higher than the transmit degree 
of freedom. By contrast, the SDP method is able to provide a secrecy rate of higher than 1.5bps/Hz even 
with K = 10. 

2) Secrecy rates versus Eves' received signal strength: We investigate the impact of 'near-far' effects 
on the secrecy rate behaviors; i.e., how the various methods perform when Eves' received signal strength, 
characterized by pi, changes. Fig. |2tb) shows the secrecy rates of the various methods with respect to 
Pg. Note that pe > 1 means that every Eve has a stronger received signal strength than Bob, while 
Pe < 1 means the vice versa. The following two phenomenons are observed for the MRT methods: First, 
plain-MRT approaches SDP for small p^, which makes sense since one may simply ignore weak Eves 
in the transmit design. Second, projected-MRT has its secrecy rate invariant to pi, which is due to its 
nulling process. This means that projected-MRT can cope with strong Eves rather effectively. Fig. [2l^b) 
also illustrates that the proposed SDP methods provide better performance than the two MRT methods, 
especially for pi > 0.5. 



10 




-10 10 20 

Transmit power (dB) 



Fig. 3. Secrecy rate versus transmit power. 

3) Secrecy rates versus the transmit power: We are interested in evaluating the secrecy rate perfor- 
mance of the various methods with respect to the transmit power budget P. The results are displayed in 
Fig. [3] Interestingly, we see that the secrecy rate of SDP appears to be approached by that of plain-MRT 
for small P, and by that of projected-MRT for large P. 
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B. The Imperfect CSI Case 

The simulation settings in the imperfect CSI case are generally identical to those of the perfect CSI 
case above, except that the average power limit is increased to P = 20dB. Regarding the imperfect CSI 
effects, we define the following channel uncertainty ratios: 



at 



A: = 1, 



and use them to control the amount of channel uncertainties in the simulations. We fix ae,i = ■ • ■ = 
ae,K — OLe- We will choose a\, = 0.03, Oe = 0.1, unless specified. The performance measure is the 
worst-case secrecy rate defined in (fT7])- ((2TI) . This performance measure does not have a closed form, but 
can be computed via SDP; the details are described in Appendix 10 

We evaluate the performance of the robust SDP method developed in Section ITV-C[ the non-robust SDP 
method (in Section IIII-BI) . projected-MRT, and plain-MRT. The latter three methods use the presumed 
CSIs h, Gi, . . . , Gk to perform transmit designs in the simulations, and then we evaluate the resultant 
worst-case secrecy rates. 



— e— SDP (proposed) 

- ► - Non-robust SDP 

- B - projected-MRT 

- A - plain-MRT 
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-S — robust SDP (proposed) 

■ > - non-robust SDP 

■ B - projected-MRT 

■ A - plain-MRT 
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(a) 



(b) 



Fig. 4. Secrecy rate performance with imperfect CSI. (a) Worst-case secrecy rate versus Eves' channel uncertainty ratio; (b) 
Probability of non-negative worst-case secrecy rate versus Eves' channel uncertainty ratio. 



1} Secrecy rate performance versus Eve's channel uncertainty ratio: Fig. |4l^a) presents the worst-case 
secrecy rates of the various methods versus Eve's channel uncertainty ratio Og- As seen in the figure, the 
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robust SDP method yields the best worst-case secrecy rate among all the methods, especially when Oe is 
large. Moreover, we observe a peculiar behavior — that non-robust SDP, projected-MRT, and plain-MRT 
yield negative worst-case secrecy rates for > 0.18. That is because these perfect-CSI-based methods 
aim to provide non-negative secrecy rate results for the presumed CSIs, but not for the actual CSIs. To 
get a better idea of how sensitive the non-robust methods can be in the presence of imperfect CSIs, in 
Fig. Ill^b) we show the probabiUty of non-negative secrecy rate; i.e., the chance that a method gives non- 
negative worst-case secrecy rate under the 1000 independent trials. One can see that the robust method 
always guarantees a non-negative secrecy rate (which is expected from its design formulation), and that 
the non-robust methods violate non-negative secrecy rate quite seriously especially for large Oe- 




Transmit power (dB) P,, 
(a) (b) 

Fig. 5. Secrecy rate performance with imperfect CSI. (a) Worst-case secrecy rate versus transmit power; (b) Worst-case secrecy 
rate versus Eves' received signal strength. 

2 } Worst-case secrecy rates versus transmit power, and Eves ' received signal strength: More results are 
shown to demonstrate the robustness of the proposed robust SDP method. Fig. [S^a) plots the worst-case 
secrecy rates of the various methods against the transmit power P. As seen, the worst-case secrecy rate 
performance of the robust SDP method is better than those of the other methods. Moreover, we observe 
that for non-robust SDP and projected-MRT, keeping on increasing P ends up decreasing the worst-case 
secrecy rate. That is because increasing the power may also help improve eavesdroppers' receptions, 
if the transmit design does not take channel uncertainties into account. Fig. [S^b) plots the worst-case 
secrecy rates against pi. Again, the robust SDP method is seen to yield better worst-case secrecy rates 
than the other methods. 
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VI. Conclusion and Discussion 

To conclude, this paper has addressed the transmit covariance design problem of MISO secrecy-rate 
maximization in the presence of multiple multi-antenna eavesdroppers, using an effective SDP approach. 
Both perfect and imperfect CSI cases are considered in our designs. We have shown by analysis that 
transmit beamforming is generally the secrecy-rate optimal strategy for the considered scenarios. As 
illustrated by the simulations, the proposed SDP solutions outperform some other existing methods. 

It is worthwhile to mention some possible extensions of this work. In B2l . we have illustrated that 
the same SDP approach can essentially be used to solve MISO SRM with more complex covariance 
constraints, such as per-antenna power constraints, and interference temperature constraints arising from 
CR. In ||43l . we have considered a multicast extension; i.e., multiple Bobs receiving the same information. 
There, our SDP approach is provably optimal only for some special cases such as the 1 Eve, 2-to-3 Bobs 
setting; it serves as a tractable approximation otherwise. Despite this, the SDP approach was empirically 
shown to provide good lower and upper bounds for the multicast secrecy rate. 

Before closing this paper, we should highlight that some concurrent research studies, e.g., |[T2l . |[T4l . 
HH, l[23l . have demonstrated that adding artificial noise (AN) in the transmit design is quite effective 
in degrading Eves' receptions. Hence, a meaningful future direction would be to extend this work to 
the AN-aided case, optimizing the transmit design and AN jointly for the maximum secrecy rate. The 
resultant design optimization is expected to be more difficult than those tackled here, and it would be 
interesting to see how the proposed SDP approach may be used to help overcome the arising design 
challenges. Moreover, while this paper has its theme on optimal SRM designs, a parallel direction that is 
also meaningful is to further investigate simple suboptimal SRM designs and their aspects. For example, 
one intuitive design is to treat all the Eves as one "super-Eve", and then employ the one-Eve optimal 
design ^ which has a closed form. Owing to page limit, such a possibility is not considered in this 
paper, but would be another interesting future direction. 
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Appendix 

A. Proof of Proposition |7] 

For Problem (ITTI) . let us first write out its Lagrangian function: 

K 

jO{W,Y,fj.) = Tr(W) - Tr(WY) + ^fik (2^(1 + ']>(GfcGf W)) - Tr(hh^W) - l) , 

k=l 

where fi = [ fii,...,fj,K ]'^, fJ-i fJ-K > are the Lagrangian dual variable for the minimum 
secrecy-rate constraints, and Y G H^* is the Lagrangian dual variables for the constraint W ^ 0. 
The corresponding KKT conditions are shown to be 

Y = I + 2^ Ef=i /^fcGfcGf - (Ef=i f^k)hh^, (35a) 

YW = 0, (35b) 

1 + Tr(hh^W) > 2^(1 + Tr(GfcGf W)), (35c) 

W^O, Y ^ 0, /Ufc > 0, k = l,...,K. (35d) 

Note that in general, problem (fTTI ) satisfies Slater's constraint qualification condition: If problem ([TTI ) 
has a feasible point, then one can prove, by construction, that there exists a strictly feasible point for 
problem (fTTI ). As a result, strong duality holds and the KKT conditions are the necessary conditions for 
a primal-dual point (W, Y,/x) to be optimal. 

The key to showing the rank-one structure of W lies in (I35al ). Let 

B = I + 2^ Y,k=i i^kGkG^. 

We see that B is positive definite, and thus has full rank. By letting p = l^k > and by denoting 

B^/^ as a positive definite square root of B, we have that 

rank(Y) = rank(B~^/2YB-i/2) 

= rank(I - p{B~^/'^h)(B-^/^h)^) > Nt - 1, 

i.e., rank(Y) is either Nt or Nt — 1. For rank(Y) = Nt, (I35bl ) can only be satisfied by W = 0. However, 
W = violates (I35cl ) when R > 0. For rank(Y) = Nt — I, (I35bl ) is achieved only when W lies in the 
nuUspace of Y, the dimension of which is one. This means that any optimal W must be of rank one. 

The proof above has shown that any primal optimal solution W of (fTTb must be of rank one for 
R > 0. Next, we consider the uniqueness of the optimal W. Suppose that there are two distinct optimal 
solutions, say Wi and W2, which satisfy rank(Wi) = rank(W2) = 1. It can be easily shown that 
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the subspaces spanned by Wi and W2 must be different in order for Wi and W2 to be distinct; i.e., 
7^(Wi) 7^ 7^(W2) where Tl{-) denotes the range space of the argument. As a basic result in convex 
optimization, any W3 = /3Wi + (1 — /3)W2, for /3 G (0, 1), is also an optimal solution ||35l . Since Wi 
and W2 are distinct, it can be easily shown that W3 is of rank two, which violates the necessity that 
any optimal W must be of rank one. In other words, we must have one optimal W only. 

B. Proof of Theorem [7] 

For ease of exposition, we restate Problems ( fT3l ) and (fT4b in the following equations, respectively 

7reZax = ™™ 0(W) 

wto (36) 

S.t. Tr(W) < P 

and 

min Tr(W) 

s.t. 7.^eZax > W 

where, with a slight abuse of notations but for notational simplicity, "f*c.iaxi-^) replaced by 7*2;^^., and 

' k=l,...,K 1 + h^Wh 

is used to denote the objective function of (l36l ). 

Our proof is divided into three steps: First, we prove that an optimal solution to Problem (|37] ) is also 
an optimal solution to Problem ( [36l l; second, we prove the converse; finally, we utilize Proposition 1 to 
establish our claim in Theorem 1. 

Step 1: Let W be an optimal solution of (l36l ). and W be an optimal solution of (|37] ). By noting that 
W is also feasible to ( [37l ). we have that 

P > Tr(W) > Tr(W). 

Hence, W is also feasible to ( [36b . This further implies that cl)(W) > 'yreiax- Moreover, as an optimal 
solution of (l37l) . W must satisfy the constraint in (l37l) ; i.e., "f*^ia^ > ipi^)- We therefore have 0(W) = 
relax''' Other words, W is optimal to ( [36l ). 
Step 2: Suppose that W is optimal to Problem ( [36l ). but not optimal to ( [37l ). Since W is feasible 
to (I37] ). the following relation holds 

Tr(W) < Tr(W) < P. 
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Since Tr(W) < P, we can construct another point W = oqW with oq > 1> such that Tr(W) = P; 
i.e., W is feasible to Problem (l36l) . Now, consider the following function 

1 + a max Tr(GfWGfc) 

f{a) := (/>(aW) = ''_zli:2± ^ . 

^ ^ ^ ^ 1 + ah^Wh 

The function f{a) is strictly decreasing with respect to a: It can be verified that 

max Tr(Gf WGfc) - h^Wh 

/'(q) = '^^^^ < 0, 

(l + ah^^Wh)2 

where the inequality above holds for h^Wh > ^ max^ Tr(G|^WGfc), which must be true for < 
7* < 1. With the strictly decreasing property of f{a) and ao > 1, we have that 

</.(W) = <A(aoW) = /(ao) < /(I) = </.(W) = 7,*,,„, 



which means that W can achieve a lower objective value than W in Problem ( [36] ). This contradicts the 
optimality of W for Problem (l36l ). 

Step 3: So far, we have proven that Problems ( [36l ) and ([37]) have the same optimal solution set. Since 
'^reiax ^ ^' Proposition 1, we know that there is a unique rank-one optimal solution to Problem ( [37] ). 
Hence Problem ([36] ) also admits a unique rank-one optimal solution. Moreover, this rank-one optimal 
solution fulfills the equality 7* = 7*g;„^. (recall Lemma 1), thereby serving as an optimal solution to the 
original problem ([T2|) as well. On the other hand, let W* be an optimal solution of ([Tit . Since W* is 
feasible to Problem ([36] ) and the relaxation is tight, i.e., 7* = 7*g;„^., W* is optimal to Problem ([36] ). 
too. This further implies that W* has to be unique and of rank one. 

C. Proof of Proposition [3] 

The key to the proof lies in the KKT conditions. The Lagrangian function of Problem ( [30] ) is given 

by 

K 



C{X) = Tr(W)- J]Tr(Tg,fc(W,Ag,fc,^)Ag,fc) 



k=l 

K 



-Tr (T,(W, \, 9) A,) - TY(WS) - X^fib - ^e,kl^e,k (38) 



k=i 



where X = |W, S, A;,, Ag, ^, /if,, /Xg, Ag, A;,} collects all the primal and dual variables, with Ag = 

[ Ag,i, . . . , Xe,K f, Me = [ /Ue,i, • • • , f^e,K f, and Ag = {Ag,,.}f^,; Ag,fc G A, G H^^'+\ S G 

fih G M_|_ and /Xg fc G M+ are dual variables associated with Tg^^, Tf,, W, A^ and Ag^^, respectively. 
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For ease of expression, we rewrite Tf,(W, A;,, 9) and Tg ^(W, Ae,fc, 9) as 

T,{W,Xb,9) = n{Xh,9) + \^WVt„ 



(39) 
(40) 



where 







-Xbel-e + 1 









-Ae,fee^^ + 2-^0-1 



Substituting (1391 ) and (1401 ) into (1381 ). we obtain an alternate expression of the Lagrangian function 

K 

C{X) = Tr(W) + Tr(WfcVe,fcAe,fcV^fc) - Tr(WVbA,Vf ) 

k=l 

-Tr(WS) + ip{Xb, Xe, 9, Hb, /^e, Ae, Ab) , 

= Tr(W) + J] J] TV(WB^'f ) - Tr(WVfeAbVf ) 

fe=i «=i 

-Tr(WS) + ip{Xb, Ae, 61, ^fe, /Ig, Ae, A;,) , 

where b|,''P G H^' is a block submatrix of Vg^fcAe^fcV^^; specifically, 



(41) 



B 



(1,1) 

e,k 



B 



(l.Afe.O 

e,fc 



B 



(Afe..,l) 



... B 



and (/?(Afe, Ae, ^, /ib, /^e; Ae, Afe) collects the terms not related to W and S, which are not important to 
the proof. 

We consider only the KKT conditions relevant to the proof here: 



VwCiX) = 0, 
Tb{W,Xb,9)Ab = 0, 
WS = 0, 

Ab > 0, W ^ 0, Abh 0, Ae,fc ^0, k = l, 



,K. 



Using the expression (HTI) . the KKT condition (I42al ) is obtained as follows: 

K N^.k 



(42a) 
(42b) 
(42c) 
(42d) 

(43) 



k=l 1=1 



April 19, 2011 



DRAFT 



27 



Premultiplying the two sides of (l43l) by W, and making use of (I42cb . we get 

K N^.k 

W (l^. + E E ^e.f ) = WV, AfeVf . (44) 

k=l 1=1 

Now the following relation holds: 



ranl< 



k(W) = rank W (l^, + ^ ^ B^'f ) (45a) 



= rank(WV6A6Vf ) (45b) 

< min{rank(V6AbVf ), rank(W)} (45c) 

where (l45ab is due to Iat^ + Ylk^i Ez=T ^fk ^ 0' <l45bl) and (l45cl) follow from (HH) and a basic rank 
inequality property BTI . If we can prove that rank(V5Af,V^) = 1, then, from ( |45l ). we will obtain 
rank(W) < 1. Therefore, in the remaining part of the proof, we will focus on the rank of V;,A;,V^. 
Substituting into the KKT condition ( |42b| ). we obtain 

Tk{Xb, e)Ah + Vf WVfeAfc = 0. (46) 

And it follows by postmultiplying (l46l ) by that 

n{Xb, AfcVf + WVfeA.Vf = 0. (47) 

By noting the following facts 

[I;v. 0]mXb,9) = Xb[lN, ] = Afe(V, - [ 0^, h]), 

[I,v. 0]Vf = In,, (48) 
we premultiply the both sides of (l47l ) by [ Itv, ] to get 

Afe(Vfe-[0^, h])AfeVf +WVfeA,Vf = 0, (49a) 
^ (AfeIjv.+W)VfeAfeVf = Afe[07v, h]AfeVf. (49b) 
We claim that A;, must be positive. Suppose that A;, = 0. Then, according to (I49al ). we have WVf,AfeV^ = 
0. By (m and In, + Ef=i Y^iti ^ek ^ 0' have W = 0. However, W = is infeasible to the 
relaxed R-SRC problem (l24b . whenever > 0. Therefore, Af, > must hold. With A;, > 0, we have that 



rank(VfeAfeVf ) = rank((AbI,v. + W)Vf,AbVf ) (50a) 
= rank(A6[ Otv, h ]AfeVf ) (50b) 
< rank([ Oat, h ]) < 1, (50c) 
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where (ISOal ) is due to A^I at^ + W ^ 0, (ISObI) and (ISOcb follow from ( |49bl ) and a basic rank inequality 
property iHTI . 

Combining (l45l) and (l50l) . we have 

rank(W) < rank(VbA;,Vf ) < 1. 

Since W 7^ for i? > 0, the rank of W must be one. 

Regarding the uniqueness of the optimal solution, the proof is exactly the same as that in Proposition 1 . 
We therefore omit it for brevity. 

D. Proof of Proposition [4] 

By the change of variable W = Z/^, ^ > 0, Problem ( [32l ) can be transformed to 

t + max max 

fc=i,...,XG.eB.,. " 

mm Tj (51a) 

z,« P + min h^Zh 

s.t. TY(Z) < iP, (51b) 
Z ^ 0, ^ > 0. (51c) 
We first show that (ISTI ) is equivalent to the following problem 

min £ + max max TrfG^^ZGt.) (52a) 

Z,5 ^ /fc=l,...,XG,eBe,fc ^ 

s.t. £ + min h-f^Zh > 1, (52b) 

Tr(Z) < eP, (52c) 
Z )r 0, £ > 0. (52d) 

Consider the optimal solution of (l52l ). say, denoted by (Z*,£*). From (I52bl ) and (I52cl ). it can be verified 
that £* > must hold. Hence, (Z*, is feasible to ([Sl])- One can deduce that (Z*, £*) is optimal to (ISTI ) 
if it holds true that 

C + min h^Z*h = 1. 

We use contradiction to verify the latter. Suppose that £* + minhg^b h^Z*h > 1. Then we can choose 
a feasible point (Z, £) = (aZ*, a£*) for some < a < 1 such that £ + minhg^t h^Zh = 1. The point 
(Z, £) can be verified to achieve an objective value lower than that offered by the optimal point (Z*, £*), 
which is a contradiction. 
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mm r 



Our next step is to turn (1521 ) to an SDR Using the epigraph reformulation, (152] ) can be rewritten as 

(53a) 
(53b) 
(53c) 

Tr(Z) < ^P, (53d) 
Z >i 0, C > 0. (53e) 
By applying the 5-procedure to convert the constraints ( I53b| ) and ( |53c| ) into LMIs, we obtain the SDP ( [34l ). 



s.t. r>^+ max Tr(G|^ZGfc), /c = 1, . . . , K, 
£ + min h^Zh > 1, 



£. Worst-case Secrecy Rate Calculation 

The worst-case secrecy rate function ipCW) in ([T8])- (|2TI ) can be computed when W is of rank one 
(which is the case of all the considered methods in this paper). With rank-one W, V'(W) can be reduced 
to 



■0(W) = minlog(l + h^Wh) - max max log(l + Tr(G^ WGfc)). (54) 

h€Bb k=l,...,K Gk&Be,k 



The first and second terms of (1541 ) are separate optimization problems. They can be recast as the following 
two SDPs by using the 5-procedure: 



Ti = max Ti 

Tl,Ab 



s.t. 



XbiN, + W 



Wh 



h^W h^Wh + 1 - n - Xbel 



Xh > 0, 



and 



r2 = mm T2 



S.t. 



h 0, 



Ae,fc>0, k = l,...,K, 

where Wfc = Iat^ ^ (g) W. Once the optimal values and T2 are computed (by an available SDP solver), 
the worst-case secrecy rate is obtained as 

V(W)=l0g(Ti*)-l0g(T2*). 
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